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ABSTRACT 



Measurements of the concentration of N 2 in the thermosphere 
obtained by the Thermosphere Probe, and temperatures derived 
from the resulting density profiles when compared with similar 
parameters derived from satellite drag, as represented by cur- 
rent model atmospheres, lead to inconsistencies. Typically, the 
near daytime maximum and the near nighttime minimum density 
observed (2. 0 x 10 9 /cc and 1.2 x 10 9 cc at 200km) and the cor- 
responding temperatures (825°K and 725° K) indicate that the 
assumption, for model atmosphere purposes, of constant boundary 
conditions at 120km should be expanded to reflect, probably, a 
significant diurnal variation. 
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FOR THE 150 TO 300 KM REGION 
AND THEIR IMPLICATIONS 
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INTRODUCTION 

Study of the orbit decay rates of sate Jlittd has led to a large amount of e- 
vealing atmospheric structure data which, on the basis of acceptable theoretical 
considerations, have been used to describe the density, temperature, composi- 
tion, and other properties of the Earth':, neutral particle atmosphere. These 
important results have appeared in various publications and also have provided 
constraints and boundary conditions for a number of model atmospheres, of 
which those of Harris and Priester, 1964, and Jacchia, 1964, are probably the 
most prominent. Also, investigators using rockets and satellites, and employ- 
ing "direct" measurement techniques with density gages, which "count" particles 
under essentially known physical conditions of volume, temperature, collection 
geometry, etc. , at rather specific points in the atmosphere, have provided a 
much smaller, but an increasing quantity of new data which permit an independ- 
ent determination of the same atmospheric parameters (Spencer, et al, 1965; 
Nier, et al, 1964; Newton, et al, 1965). Initial comparisons of these data re- 
vealed that density values derived from the drag measurements are higher than 
density values determined from the gage techniques , by approximately a factor 
of two. Temperature profiles obtained from the "direct" measurements when 
compared with drag-derived values also exhibit differences which were not con- 
sidered significant. However, a recent re-evaluation of the earlier density and 
temperature data and, in addition, consideration of new data presented in this 
paper which tend to reinforce earlier tentative conclusions, lead us to conclude 
that there are significant differences. 

We feel that the differences observed have a clear consistency that indicates 
a systematic difference exists between the results obtained by the two different 
techniques , and thus suggest either that there are appreciable (unknown) errors 
in the direct measurement data or that there exists a need for modification or 
expansion of certain of the assumptions upon which the various model atmos- 
pheres, which reflect, largely, the "drag" data, are based. Accordingly, this 
paper will present the data upon which these comments are based, an evaluation 
of the errors sustained by the gage techniques , and some supporting discussion 
of the conclusions drawn. 
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EXPERIMENTAL TECHNIQUE 

It is well known among vacuum laboratory workers that the measurement 
of the absolute concentration of a particular gas in an evacuated device to a 
precision of ‘IS-SO percent poses difficult problems particularly when the total 
pressure is in the 10" 6 torr region or below. The task is least difficult for 
chemically inactive gases such as molecular nitrogen, for which reasonably 
well -understood and reliable procedures exist- Fortunately , the most plentiful 
constituent of the lower thermosphere is molecular nitrogen This fact, and its 
chemical stability, thus make it likely that measurements oi N 2 concentration, 
using laboratory related devices , will lead to concentration data for this atmos- 
pheric constituent in which one can have high confidence . This belief is rein- 
forced by the fact that a relatively satisfactory preflight calibration can be 
accomplished for this gas, in comparison, for example, to atomic or molecular 
oxygen, for which only very questionable caljHi Jy ; ? c&n be performed with the 
present state-of-the-art techniques. 

The data to be discussed in this paper, molecular nitrogen cot oentration 
and temperature , have been obtained through application of a small c legatron 
mass spectrometer in an experimental approach which is based !&:>• J 7 upon 
kinetic theory considerations. The nerimantal technique was disc/* ied in 
detail in an earlier paper (Spencer al, 1965) and thus will only r* omnia- 
rized ve^y briefly here. 

An atmospheric composition meao orement system (see Fig> * • which has 
been called the Thermosphere Probe (TP) <•? arranged to be a. • -a aloft and 
then ejected into the thermosphere from its la>ujj.h'ns; vchic ,ure 2). The 
TP, which is tumbled intentionally end-over-end after f-.jee; • u room the rocket, 
contains the omegatron, ail supporting electronic systems, nod an appropriate 
telemetry system for N 2 data transmissions to ground si&tio.ts. Data trans- 
mitted from the instrument also includes information relative to its orientation 
which is essential for satisfactory data analysis, and various additional data 
such as sensor temperature. A representative telemetry record showing N 2 
pressure variations occurring in the omegatron chamber as the TP tumbles 
through the atmosphere is shown in Figure 3. 


RESULTS OF PREVIOUS EXPERIMENTS 

The data analyses procedures employed to derive the desired atmospheric 
Nj density and temperature profiles from the observed omegatron chamber 
pressure, and other data, have been discussed in detail in the earlier publica- 
tion. In brief review, the density profile is directly determined from measured 
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Figure 1 . Schematic Representation of Thermosphere Probe (TP) 



Figure 2. Artiit's Representation of Thermosphere Probe Ejection Procedure 


3 




ALTITUDE (Km) 





4 


Figure 3. Photograph of Telemetry Record of Typical Omegatron Chamber Pressure 



ALTITUDE (kilometers) 


gage density and instrument orientation and velocity. The temperature data are 
then derive^ from the density profile, employing scale height considerations. 
Atmospheric data resulting from these computations are shown in Figures 4 and 
5. Harris and Priester 1964 model values most closely approximating solar 
activity at the times of the experiments are also shown to provide a basis for 
comparison. A summary of flight statistics of all launchings conducted to date 
appears in Table I. 



Nj NUMBER DENSITY (portlcln/cc) 


Figure 4. Summary of Measured N z Density Profiles Shown with 
Appropriate Model Atmosphere Profiles 
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Figure 5. Summary of Early N 2 Temperature Profiles with 
Appropriate Model Atmosphere Profiles 
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There are a number of comments pertinent to these density and temperature 
data: 

1. All profiles but those labelled NASA 6.08 pertain to the "normal" atmos- 
phere. The 6.08 data were obtained during the July 20, 1983 s( lar 
eclipse and will be discussed separately. 

2. The general similarity in shape of both the N 2 density and temperature 
data to model profiles confirm gross hydrostatic equilibrium throughout 
toe altitude region. The values of temperature, and the altitudes at 
which isothermy is observed are in general agreement with accepted 
concepts; however, the data generally indicate isothermy to significantly 
lower altitudes than the model, suggesting possibly that there are fac- 
tors ether than conductivity or mass transport that need consideration 
in model studies. On the other hand, more recent data shown and dis- 
cussed later in the paper support the gradients depicted by the model. 

3. The density data from these earlier flights which are believed to have an 
absolute accuracy of about 50% or better show a systematically lower 
density (N 2 ) atmosphere than depicted by the model. In addition, the 
density values at 150 kms, corresponding to various conditions, show a 
variation that is inconsistent with the common assumption (for model 
purposes at least) that the atmosphere at 120 km is essentially invarient. 

4. The temperature of the atmosphere is somewhat lower than indicated by 
the model as illustrated by the NASA 6.07 temperature profile. There 
are a number of considerations bearing upon this conclusion which will 
be discussed later in the paper. 

5. The profile (6.09) corresponds to winter night, quiet conditions and 
should reflect nearly minimum solar cycle density for the time (2209 
local time) and latitude (Wallops Island 38° N) , although the measure- 
ment was carried out some six hours prior to the nominal 0400 minimum 
time. We did have concern, however, over the seemingly low absolute 
density observed although there are no known aspects of the flight or 
experiment which lessen our confidence in the data. To confirm the re- 
sult however, a later experiment was arranged, as will be discussed 
below. 

6. The density of the atmosphere during the solar eclipse of July 20, 1963 
was notably different than "normal, " as reflected in the 6. 08 profile. 

The 6.07 density profile provides a convenient reference as it was ob- 
tained at about the same local time and level of solar activity three 
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months earlier. The chief difference in conditions was that a p for the 
6. 07 experiment was about 30 as compared with a value of about six 
during the eclipse. One should expect therefore, that the 6.07 density 
should be greater than 6. 08 in the 150km region, but this was not ob- 
served, and the reason may lie in the absolute error of the measure- 
ments. It does not appear to be related to semi-annual or seasonal 
variations, which would favor a lower rather than higher density. Ig- 
noring the difference at the 150 km level (absolute value difference) , and 
noting that the relative accuracy of the data is much better (probably 10% 
or less) than their absolute accuracy because absolute calibration is not 
a significant factor, the density at 300km decreased substantially during 
the eclipse. This is borne out by the 6.08 temperature profile, which 
resembles a nearly nighttime exospheric temperature condition. It is 
important to note, however, that the temperature data derived using 
scale height considerations depend upon the maintenance of an adequate 
state of diffusive equilibrium in the atmosphere , and thus the eclipse tem- 
peratures maybe less accurate than, for example, the 6.07 temperatures. 

Assuming, however, that the data are essentially correct, the low tem- 
perature reached by the time of eclipse , and thv, isothermy observed 
down to about 200 maximum km, reflect the substantial decrease of solar 
UV and support conclusions of Nicolet (i960) rega^di g conduction and 
conduction times. The fact that the density at this time is significantly 
lower than that considered by Nicolet (nearly a factor of five) indicates that 
the conclusions he reaches probably apply, at this time in the solar cycle, 
to significantly lower altitudes (1-2 scale heights) . This may explain why 
the temperatures during the eclipse decreased to values which were to us, 
surprisingly low. 


NEW EXPERIMENTS 

The experiments and data discussed above have shown densities and ten • 
peratures systematically lower than the model values chosen for comparison. 
They have also revealed other aspects , such as relatively large variations of 
gradients in the 150 - 200km region. To help resolve this and other questions, 
several steps have been taken: (a) a pair of experiments to provide a daytime 
near-maximum N 2 density measurement, and a nighttime near-minimum N 2 
density measurement, have been carried out in a manner to assure nearly iden- 
tical conditions and with identical instruments, simultaneously calibrated; (b) 
a re-evaluation of the errors encountered in the experiment has been conducted; 
and (c) a review of the analysis of the earlier experiment data discussed above 
has been initiated to reaffirm or reveal the errors sustained. 
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The results from the two flight experiments are now available and will be 
discussed, and a preliminary discussion of the error evaluation appears in the 
appendix. Re-evaluation of the analysis of the earlier data is still underway 
and will be reported in a subsequent paper . 


DISCUSSION OF NEW DATA: DENSITY 


The density data obtained from the day-night pair of experiments is shown 
in Figure 6, with model values appropriate to existing conditions, linearly 
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extrapolated (the logarithm) to the time of the experiment. The two profiles 
represent very nearly the maximum and minimum N 2 density under conditions 
of essentially quiet sun (10.7cm flux ~75xl0 -22 w/m 2 c/s). The absolute 
accuracy of the concentrations ie believed to be 25% or better and the relative 
accuracy over the altitude range, and between profiles is believed to be 10% or 
better. The following comments refer to these data: 

1. The data show a lower density N? atmosphere at 250 km than depicted by 
the model by a factor of two, with the difference gradually decreasing 
with decreasing altitude as expected. A reasonable extrapolation of the 
data for a few km to 150 km shows that the data are lower than the model 
values by about 40% of the model value at this altitude. 

2. The ratio of daytime -to-nighttime density is significantly different than 
that of the model values . These ratios (model ratio corrected to flight 
local times) are shown in Figure 7, where, it should be noted, the data 
ratio curve has been "eyeball"-extrapolated below 170 km. The large 
difference in ratio seen in the 150-200 km region is of particular interest 
and will be discussed below. 

3. Nier's data (Nier, et al, 1964), although somewhat below the nighttime 
profile at the lower altitudes, show higher values near 200km as should 
be expected (dawn vs nighttime) . Conditions at the time of the measure- 
ments were made were: time - 0730, location - White Sands, F = 77, 
a p * 16. There is good agreement with the data presented here. 


DISCUSSION OF NEW DATA: TEMPERATURE 

The temperature profiles of the thermosphere derived from the densities 
shown in Figure 6 are presented in Figure 8, with model atmosphere profiles 
appropriate to the conditions at the time of the experiments (Harris and 
Priester, 1964). The model values have also been extrapolated to the local 
times of the launchings which, as has been noted, correspond very nearly to 
the times of expected maximum and minimum temperature. Also plotted for 
comparison are two temperature curves of Jacchia (1964) , selected on the basis 
of exospheric temperatures most closely in agreement with the measured values. 
Error flags are shown on the derived N 2 temperatures , which require a brief 
discussion. It was stated earlier that the density profiles were believed accu- 
rate (relative) to 10% and on an absolute basis to 25%. Temperature calculated 
using scale height as the basis depends only upon the relative accuracy (the 
"shape") of the density profiles. The error flags on the temperature data, 
therefore, reflect the 10% flag on density, with, howevez 1 , the reservation that 
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not all possible curves that could be drawn within the density error flags are 
acceptable. Clearly, only smoothly varying ■ rofiles are acceptable, such as 
through the maxima, the minima, the mean, or other reasonably acceptable 
points. In this case, the extremes were determined on the basis of density 
profiles drawn through the high altitude maximum and low altitude minimum 
value and conversely, each passing through the medium altitude mean density. 
These choices should provide extreme temperatvre profiles. The flags on the 
temperature curves thus correspond to density curves that reflect the maxima 
and minima values appropriate to the stated relative accuracy of the density . 

There are difficulties, chr’lenges, or dangers in providing a meaningful 
discussion of the measured temperatures, however, it is probably most fruitful 
to discuss the temperature profiles in comparison with the model values of 
Harris and Priester, which for comparison here were chosen on the basis of 
known solar activity, and the model values of Jacehia which for comparison 
here were chosen on the basis of measured exospheric temperature . One of the 
reasons that makes comparisons somewhat uncertain is that the measured 
values have probable errors which at least permit consideration of improbable 
gradients. Another reason is that the models must necessarily be considered 
limited, because not only of probably over-simplified assumptions, but because 
of tne lack of realistic boundary conditions at 120km. The measured tempera- 
tures are, for example, derived from densities which show a diurnal variation 
of a factor of 1.7 at 200km which is not compatible with either model. Harris 
and Priester, for example, as shown in Figure 7, have a ratio of 1.2 at 200km. 

In spite of these difficulties , we feel it is useful to note the apparent differ- 
ences , as it may lead to a better understanding of the pertinent processes and 
may contribute to bettering the model representation of the real atmosphere. 

J. The measured exospheric diurnal temperature ratio is about 1. 15 in 
contrast to the value 1 . 20 (determined for time of day and year the ex- 
periments were conducted, from Figure 6a, SAO Special Report #150, 
April 1964, Jacchia) established by a considerable amount of satellite 
drag data. It must be noted, also, that the error flags permit temper- 
atures which have a ratio of 1.20 or greater. Also, the daytime maxi- 
mum (measured) is much less than the Harris and Priester model value, 
although showing isothermy to a comparable altitude. The Jacchia 
"850° K" curve, on the other hand, is isothermal to an altitude lower by 
50 km, as might be expected (less energy absorbed above 250 km). 

2. The measured exospheric tempera' ^re values are in better agreement 
with Jacchia temperature data (Figure 8) as deduced from satellite drag 
using, apparently, a revised Nicolet model. It is interesting to compare 
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the temperatures with those deduced by Nicolet (1963) for minimum 
solar activity conditions which occurred in June 1956. See Table II. 



Day 

Night 

Experimental 

825 

725 

Nicoiet/1963 

836 

629 


.able II 


3. Although the measured nighttime exospheric temperature agrees reason- 
ably well with the model values (Figure 8) , the gradients are not in as 
good agreement, particularly at the 250km level. In contrast to the 
nighttime profile, however, is the 6.09 result discussed earlier in the 
paper, and shown in Figure 5, where isothermy is observed at a consid- 
erably lower altitude. No solar or magnetic activity effects are known 
to account for the difference, and no explanation for the difference is 
apparent at this time, unless it is due to variable and unknown conditions 
below 120 km. 

4. In considering the measured day-night profiles , one concludes that there 
is an incompatibility in respect to gradients. For example, the nighttime 
gradient at some altitudes is quite high in comparison with the corres 
ponding daytime gradient. It should be noted, however, that the insensi- 
tivity of the density profile to changes in temperature in this region 
(hence, the "large" temperature error flags) works against the estab- 
lishment of meaningful gradient detail. The magnitudes of the two pro- 
files, however, are significant and can be supported . Since a large 
diurnal density change is observed at 200 km, in which we have consider- 
able confidence and it is not accompanied by a substantial increase in the 
temperature at 200 km (Harris and Priester model requires over 200* K 
for only a 20% density increase) some other effect must be determined. 

It seems to us, therefore, that the best explanation is an increase in the 
density (and/or temperature) above the mesopause (near or below 120 km), 
an effect not incorporated in the models. Mass transport horizontally 
may also be a possible explanation, although Harris and Priester (1964) 
h ive tentatively concluded that it is not significant. 
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SUMMARY AND CONCLUSIONS 


New measurements of the N 2 density of the thermosphere corresponding to 
nearly maximum and minimum of the diurnal variation, have been obtained under 
quiet sun and low magnetic activity conditions. The density datr. have been in- 
terpreted in terms of the neutral particle temperature. 

The results permit the following conclusions: 

1. The density variation observed near 300 km is equivalent to that depicted 
by the Harris and Priester 2964 model (—2.5) but greatly exceeds that 
indicated by the model at 200 km. 

2. The diurnal exospheric temperature variation is considerably less tha:'_ 

that indicated by the same model. } 

3. These two factors (1), (2) are not consistent with present model consid- 
erations, and in particular, probably require reconsideration of the 
constant boundary values at 120 km, as both a density and temperature 
increase in the region below and near 120km could explain the data. 

4. The measured exospheric temperatures (ignoring error flags) generally 
are in agreement with the satellite -drag determined temperatures of 
Jacchia. 

5 . Improved measurement accuracy is required to establish meaningful 
details of the temperature gradient in the 200km region. 

6. Measured N 2 densities in the thermosphere are systematically lower 
than satellite-drag determined densities. A factor of two is observed at 
250 km. 
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Appendix 


ERROR DISCUSSION 


The data presented have been reduced from telemetry records similar to 
that shown in Figure 3. The deflection versus time outputs are read using 
standard techniques and are reduced to gage number density, wall temperature, 
etc . , using pre-flight calibration of the various instruments . The ambient 
number density is obtained using the equation: 


n 


a 


An. U. 

1 1 

2 ft V cos a 


( 1 ) 


where 

n B = N 2 ambient number density 

An. = measured maximum-minus-minimum N 2 number density 

U. = y r 2KT~An = most probably thermal speed of particles inside 
measuring chamber 

V = speed of TP with respect to earth 

a - minimum angle between velocity vector of TP and plane of tumble. 


The basic assumptions in the derivation are (a) the measured gas has a 
Maxwellian distribution of velocities , both inside and outside the measuring 
chambers; (b) the particles inside the chamber are fully accommodated to the 
chamber wall temperature, (c) the atmosphere rotates with the earth, and (d) 
certain implicit assumptions concerning gage geometry and homogeneity which 
will be discussed in a later section. 

As can be °een from Equation (1), the uncertainty in the calculated ambient 
number density depends upon the uncertainty in the measured quantities anj , 

U| , V and a approximately in the following manner: 
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( 2 ) 


t 


Sn a 8 (An.) 


The gage sensitivity used to obtain an* is determined from a vacuum calibration 
of the sensor. In general, in the calibration several Bayard-Alpert ionization 
gages are used as secondary standards after having been calibrated against a 
McLeod gage. Several calibrations of the omegatron are performed and self- 
consistency Ijetween the BA gages and the omegatron calibration is carefully 
established. In the case of NASA 18.01 and NASA 6.11, the two omegatrons 
were calibrated simultaneously on the same system. All of the other data from 
other flights presented were also obtained from omegatrons which were cali- 
brated on the same vacuum system with the same BA gage secondary standards . 
As a result, it is believed that the relative accuracy of the calibration data is 
better than ±10 per cent. For NASA 18.01 and NASA 6.11 the relative accuracy 
of the calibration data is probably better than ±5 per cent. An evaluation of the 
factors knov’i to effect the validity of the absolute calibration leads to the con- 
clusion that the absolute accuracy is ±25 per cent. 

Uj, the most probable thermal speed of the particles within the measuring 
chamber, is calculated from the measured chamber wall temperature, which is 
measured by a thermistor and is believed known to be better than ±2 per cent. 
Since Uj involves the square root of the temperature, the resulting error con- 
tribution in Equation (2) is on the order of ±1 per cent. 

The TP velocity, V; angle of attack, a; and the basic independent variable, 
altitude, are determined using trajectory data obtained from radar. These data 
were obtained, at Wallops Island, Virginia, from the FPS-16 and FPQ-6 radars. 
The FPS-16 specification accuracy is ±5 yards in range and ±0. 1 mil in azimuth 
and elevation. The FPQ-6 specification accuracy is ±5 yards in range and 
±0. 05 mil in azimuth and elevation. A theoretical trajectory is fitted to the 
radar data, and considered adequate when the downleg low altitude information 
agrees better than ±0 . 5 km. The resulting velocity data are better than ±1 per 
cent. Since the TP and rocket motor separate at approximately 4 feet per sec- 
ond and it is uncertain which portion the radar is tracking, ±0.5 km is adopted 
as the uncertainty in altitude at approximately 120 km altitude on the downleg 
portion of the trajectory. This is equivalent to approximately ±0.25 km at peak 
altitudes and a very small error for upleg data. This uncertainty in altitude 
implies a ±5 per cent maximum uncertainty in the density versus altitude profile 
for the low altitude downleg data. 

The angle of attack data are obtained through an optical technique using the 
method described by Taeusch, Carignan, Niemann and Nagy (1965). This 
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method requires two measured inputs to determine an angular momentum vector. 
The angular momentum vector and the velocity vector (from radar data) -?.re then 
used to determine the angle of attack. The total effect of errors in the three 
measured quantities varies greatly with geometry and is best considered for 
eajh individual case . A worst case analysis for the data presented herein yields 
a possible error in a of approximately 5° . The effect of an uncertainty in a on 
the computed density, depends on a (Equation 2). Figure 9 shows this depend- 
ency. For a stable rocket, a minimum is near zero at ejection, and typically 
less than 20° during most of the ascending leg. The ingle of attack versus alti- 
tude for NASA 18.01 is plotted in Figure 10 for illustrative purposes. For an 
uncertainty in a of ±5° at an a of 20° , the resultant uncertainty in density is ±3 
per cent. A plot of the effect of SaonSn 0 Ai a is shown in Figure 9. 
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Figure 10 . Error in Ambient Density Due to Error in a Minimum, for Various oc 






Telemetry and reading errors which contribute to the uncertainty in An, 
and Uj are on the order of ±2 per cent. Therefore, we see in Equation (2) that, 
for relative considerations: 


S(Ani) 


An. 


±10% relative calibration 
~ ±2% telemetry and reading 
±25% absolute calibration 


SU ±1% absolute calibration 

1 

u - 

* ±2% telemetry and reading 


-y- - ±1% 


or 


tan a8a ~ ±3% 


®W 


— - ±10.2% RMS relative - ±25.1% RMS absolute 


and 


SU. 


y ~ ±2.2% RMS absolute 


5n . 

~ ±10.9% RMS relative ~ ±25.4% RMS absolute 


Finally, it is appropriate to discuss briefly a correction that is made nec- 
essary by imperfections in gage geometry. As stated previously, the derivation 
of Equation (1) required certain assumptions concerning gage geometry and 
"homogeneity" of the particles within the gage. One major assumption made in 
the use of Equation (1) is that the orifice diameter is small compared to the 
mean free path and the chamber dimensions. Since, however, this would, if 
strictly true, result in poor conductance between the chamber and atmosphere, 
a chamber pressure time response problem would exist . Rather than create 
this problem for the measurement, the gage orifice is made large enough to 
allow quick response to the density changes due to tumble and altitude changes. 
The output varies from that theoretically expected during a tumble cycle. 
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Therefore , empirical correction factors are determined from the data during 
those times of low minimum a, as shown in Figure 11. These corrections are 
then applied to the data obtained when the minimum a is high. Figure 12 shows 
the correction factors for various s 0 ’s for flight NASA 18.01. The effect of this 
correction factor on the data uncertainty is a function of the minimum angle of 
attack and the ability to determine the factor adequately from flight data. It has 
been determined that the maximum correction to be made to the data is of the 
order of 18%, and the correction is determinable from flight data. A negligible 
relative error is believed and assumed to result, and accordingly, is not con- 
sidered in the error calculation. 

The normalization of the measured and theoretical curves at a = 0 implies 
that the relationship expressed by Equation (1; is correct at that point. The 
validity of this assumption is now being studied and the results will be reported 
at a later time. 
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Figure 1 1 . Theoretical and Experimental Variation of Pressure Due to f(s) 
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